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' Because of two 3a anomalies, the Standard Model (SM) fit of the precision elec- 

Q^' troweak data has a poor confidence level, CL = 0.010. Since both anomalies involve 

p • challenging systematic issues, it might appear that the SM could still be valid if the 

anomalies resulted from underestimated systematic error. Indeed the CL of the global 
fit could then increase to 0.65, but that fit predicts a small Higgs boson mass, mn = 43 
■ GeV, that is only consistent at CL = 0.035 with the lower limit, mn > 114 GeV, es- 

}—i ' tablished by direct searches. The data then favor new physics whether the anomalous 

■ ■ ■ ' measurements are excluded from the fit or not, and the Higgs boson mass cannot be 

predicted until the new physics is understood. Some measure of statistical fluctua- 
tion would be needed to maintain the validity of the SM, which is unlikely by broad 
statistical measures. New physics is favored, but the SM is not definitively excluded. 
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1. Introduction 

A decade of beautiful experiments at CERN, Fermilab, and SLAC have provided in- 
creasingly precise tests of the Standard Model (SM) of elementary particle physics. The 
data is important for two reasons: it confirms the SM at the level of virtual quantum effects 
and it probes the mass scale of the Higgs boson, needed to complete the model and provide 
the mechanism of mass generation. In the usual interpretation the data is thought to con- 
strain the Higgs boson mass, mn, most recently with mn < 195 GeV[|I| at confidence level 
CL = 95%. At the same time direct searches for the Higgs boson at LEP II have established 
a 95% CL lower limit, mn > 114 GeV.gQ 

Recently the agreement of the precision data with the SM has moved from excellent to 
poor. For the global fits enumerated below, the confidence level has evolved from 0.45 in the 
Summer of 1998,0 to 0.04 in the Spring of 2001,]^ and then to 0.010 in the current Spring 
2002 data. The current low CL is a consequence of two 3a anomalies, together with the 
evolution of the W boson mass measurement, as shown below. The 3a anomalies are (1) 
the discrepancy between the SM determination of x\y = sin^9\y, the effective leptonic weak 
interaction mixing angle, from three hadronic asymmetry measurements, x\y [Ah] , versus its 
determination from three leptonic measurements, a;^y[y4i], and (2) the NuTeV measurement 
of charged and neutral current (anti)neutrino-nucleon scattering, [|| quoted as an effective 
on-shell weak interaction mixing angle, a;^y^[*'z/^ A^]. 

If either anomaly is genuine, it indicates new physics, the SM fit is invalidated, and 
we cannot use the precision data to constrain the Higgs boson mass until the new physics 
is understood. However both anomalous measurements involve subtle systematic issues, 
concerning experimental technique and, especially, nontrivial QCD-based models. If the 
systematic uncertainties were much larger than current estimates, the CL of the global fit 
could increase to as much as 0.65, as shown below. It is then possible to imagine that the 
SM might still provide a valid description of the data and a useful constraint on the Higgs 
boson mass. 

We will see however that this possibility is unlikely, because of a contradiction that 

emerges between the resulting global fit and the 95% CL lower limit, itlh > 114 GeV. The 

central point is that the anomalous measurements are the only mn -sensitive observables that 

N.B., the experimental 95% lower limit from the direct searches does not imply a 5% chance that the 
Higgs boson is lighter than 114 GeV; rather it means that if the mass were actually 114 GeV there would 
be a 5% chance for it to have escaped detection. The likelihood for niH < 114 GeV from the direct searches 
is much smaller than 5%. See for instance the discussion in section 5 of 

^ CL = 0.010 for Spring 2002 is from a fit specified below that uses the same set of measurements as were 
included in the quoted 1998 and 2001 fits. Reference ||l| has a slightly different value for their all-data fit, 
CL = 0.017, because of two recently introduced measurements, which we do not include as discussed below. 
Furthermore, updating the all-data fit of ||l| we find, as discussed below, that it would now yield CL = 0.009. 
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place the Higgs boson mass in the region allowed by the searches. All other mji^-sensitive 
observables predict mn far below 114 GeV. We find that if the anomalous measurements are 
excluded, the confidence level for mH > 114 GeV from the global fit is between 0.030 and 
0.035, depending on the method of estimation. 

The hypothesis that the anomalies result from systematic error then also favors new 
physics, in particular, new physics that would raise the prediction for nin into the exper- 
imentally allowed region. This can be accomplished, for example, by new physics whose 
dominant effect on the low energy data is on the W and Z vacuum polarizations (i.e., 
"oblique" [0]), as shown explicitly below. Essentially any value of mn is allowed in these fits. 

It should be clear that our focus on the possibility of underestimated systematic error 
is not based on the belief that it is the most likely explanation of the data. In fact, the 
situation is truly puzzling, and there is no decisive reason to prefer systematic error over 
new physics as the explanation of either anomaly. Rather we have considered the systematic 
error hypothesis in order to understand its implications, finding that it also points to new 
physics. 

The SM is then disfavored whether the experimental anomalies are genuine or not. The 
viability of the SM fit and the associated constraint on tuh can only be maintained by 
invoking some measure of statistical fluctuation, perhaps in combination with a measure 
of increased systematic uncertainty. This is a priori unlikely by broad statistical measures 
discussed below, but it is not impossible. The conclusion is that the SM is disfavored but not 
definitively excluded. A major consequence is that it is important to search for the Higgs 
sector over the full range allowed by unitarity, ^ as, fortunately, we will be able to do at the 
LHC operating at its design luminosity. 

This paper extends and updates a previously published report, based on the Spring 
2001 data set, which focused exclusively on the m^-sensitive observables. The present anal- 
ysis is based on the Spring 2002 data, and considers m/^-sensitive observables as well as 
global fits of all Z-pole observables. The data has also changed in some respects: the 3a 
NuTeV anomaly is a new development and the discrepancy between the hadronic and lep- 
tonic determinations of x^/ has diminished from 3.6 to 3. Oct. However the other m^-sensitive 
observables are unchanged, and the present conclusions are consistent with the previous re- 
port. 

Since in this work we also consider global fits, we can summarize the conclusion quan- 
titatively by introducing the combined probability 

Pc = CL(Global Fit) x CL{mH > 114). (1.1) 

The internal consistency of the global fit and its consistency with the search limit are inde- 
pendent constraints, so the combined likelihood to satisfy both is given by Pc- We find that 
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Pc is roughly independent of whether the three hadronic front-back asymmetry measure- 
ments are included in the fit, although the two factors on the right hand side of eq. (1.1) 
vary considerably in the two cases. For instance, for the global fit to 'all' data, we have 
CL(Global Fit) = 0.010 and CL{mH > 114) = 0.30 so that Pc = 0.010 x 0.30 = 0.0030. If 
the three hadronic asymmetry measurements are omitted we have instead CL(Global Fit) = 
0.066, CL{mH > 114) = 0.047, and Pc = 0.066 x 0.047 = 0.0031. The extent of the agree- 
ment in this example is accidental, but the point remains approximately valid: if the three 
hadronic asymmetry measurements are omitted, the increase in the global fit confidence level 
is approximately compensated by a corresponding decrease in the confidence level that the 
fit is consistent with the direct search limit. 

In section 2 we review the data used in the fits, with a discussion of how it has evolved 
during the past few years which emphasizes the importance of the W mass measurement. 
In section 3 we briefiy discuss the three generic explanations — statistics/systematics/new 
physics — of the discrepancies in the global SM fit. In section 4 we review the methodology 
of the SM fits and the choice of observables. In section 5 we present fits of the data which 
exhibit the range in rriH preferred by the m//-sensitive observables, as well as global fits with 
and without the anomalous measurements. In these fits we use the distribution for the 
global fits and the Ax^ method to assess the consistency of the fits with the direct search 
lower limit on mn- In section 6 we use a "Bayesian" maximimum likelihood method instead 
of Ax^ to estimate the CL for consistency with the direct searches. Section 7 illustrates 
the possible effect of new physics in the oblique approximation. The results are discussed in 
section 8. 

2. The Data 

We consider 13 Z-pole observables and in addition the directly measured values of 
mw, the W boson mass, m^, the top quark mass, Aa^, the hadronic contribution to the 
renormalization of the electromagnetic coupling at the Z pole, and the NuTeV result. As 
discussed in section 4, we do not include the W boson width or the Cesium atomic parity 
violation measurement, which is the principal reason for the small differences between the 
global fits presented here and in |^. These measurements have only recently been added to 
the global fits; they were not included in the 1998 and 2001 fits|§|, H which we also consider 
below. Our all-data fit is tabulated in table 2.1, with the current preliminary experimental 
values from [|l|. Details of the fitting procedure are given in section 4. 

The central value for x^[''z/^ A^] from the NuTeV experiment is shown in table 2.1. In 
our SM fits we include the small dependence of x^[''z/^ A^] on and rriH given in ||^. Table 
2.1 also contains the model independent NuTeV result @], given in terms of effective Zqq 
couplings, gj^ = glj^ + g'j^j^ and g\ = glj^ + g^^. They are not included in the SM fits but are 
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used instead of Xy^ [v N] m. the new physics fits of section 7. 

The confidence level of the SM fit in table 2.1 is poor, CL = 0.010, with x^/N = 27.7/13. 
The central value of the Higgs boson mass is mn = 94 GeV. As shown in section 4, our results 
agree very well with those of [|l| when we fit the same set of observables. 

The SM fit was excellent in 1998 and has now become poor. Large discrepancies occur 
among the six SM determinations of the effective leptonic weak interaction mixing angle, 
x^^. The three leptonic measurements, Alr, Ap^, and Ae^r are quite consistent with one 
another. They combine with /N = 1.6/2, CL = 0.45, to yield 



The three hadronic measurements are also mutually consistent and combine, with x^ /N = 
0.03/2 and CL = 0.985, to yield 



But 

^vi/[^i'] ^w[Ah] differ by 2.99cr corresponding to CL = 0.0028. Combining (2.1) and 
(2.2), the result for all six measurements is x\y = 0.23149(17). The very small associated 
with the three hadronic measurements is either a fluctuation or it suggests that the errors 
are overestimated, in which case the discrepancy between xjyfA^] and would be even 

greater. 

The discrepancy between x\y [Al] and x\y [Ah] is driven by the difference of the two most 
precise measurements, Alr and A'^p^, which has been a feature of the data since the earliest 
days of LEP and SLC. At present, a;|y from A^r and A^p^ are respectively 0.23098(26) and 
0.23218(31). They differ by 2.97ct, CL = 0.0030, and combine to yield x\y = 0.23151(20). 

Combining all 6 measurements directly we flnd x\y = 0.23149(17) as above, with x^ = 
10.6/5 and CL = 0.06. Notice that the ratio of this confldence level to the confldence 
level, CL = 0.003, for versus x^yfAj^], 0.06/0.003 = 20, is just the number of ways 

that two sets of three can be formed from a collection of 6 objects. If one attaches an 
a priori signiflcance to the leptonic and hadronic subsets, then the appropriate confldence 
level is 0.003, from the combination of a;|4/[74i] and ^(^/[Aj^]. If instead one regards the 
grouping into xjyfAi] and a;'^y[Aj:/] as one of 20 random choices, then 0.06 is the appropriate 
characterization of the consistency of the data.0 In either case the consistency is problematic. 

The determination of x\y from the hadronic asymmetries assumes that the hadronic 
Zqq interaction vertices are given by the SM. For instance, to obtain x\y from 



x'^[Al] = 0.23113(21). 



(2.1) 



x^^[Ah] = 0.23220(29). 



(2.2) 



A'' 



FB — 




(2.3) 



'^I thank M. Grunewald for a discussion. 
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we assume that is at its SM value, Af, = ^^[SM]. y4b[SM] has very httle sensitivity to the 
unknown value of mn, and not much sensitivity to the other SM parameters either. x\y is 
then obtained from A^ = {qIl - geR)/{9eL + 91r), using ^i^^ = ^ + and QeR = x\^. The 
only assumption in obtaining xjy from the leptonic asymmetries is lepton flavor universality. 

The 3a discrepancy between x'^^[y4i] and is significant for three reasons. First, 

it is a failed test for the SM, since it implies Ag ^ Ag[SM]. For instance, A^ extracted from 
A'^p^ (taking Ai from the three leptonic asymmetry measurements) disagrees with y4b[SM] 
by 2.9(7, CL = 0.004. Second, together with the niw measurement, the ^(^/[Al] - 
discrepancy marginalizes the global SM fit, even without the NuTeV result. Finally, in 
addition to the effect on the global fit, it is problematic that the determination of the Higgs 
boson mass is dominated by the low probability combination of ^(^[Ai,] and or 
by the low probablility combination of the six asymmetry measurements. In judging the 
reliability of the prediction for rriH we are concerned not only with the quality of the global 
fit but also with the consistency of the smaller set of measurements that dominate the rriH 
prediction. 

To understand the effect on the global fit it is useful to consider the evolution of the 
data from 19980 to the present[||, shown in table 2.2, together with the intervening Spring 
'01 data set0, on which was based. The - discrepancy evolved from 

2.4a in '98 to 3.6a in Spring '01 to 3.0a in Spring '02. Excluding NuTeV, the CL of the set 
of measurements listed in table 2.1 evolved during that time from a robust 0.46 to 0.04 to 
0.10.0 

The decrease in the global CL is only partially due to the changes in the asymmetry 
measurements. An equally important factor is the evolution of m-\y, for which the precision 
improved dramatically, by a factor of 3, while the central value increased by |a with respect 
to the '98 measurement. To understand the role of m^, table 2.3 shows fits based on the 
current data plus two hypothetical scenarios in which all measurements are kept at their 
Spring '02 values except mw- In the first of these, mw is held at its '98 central value and 
precision. In the second the current precision is assumed but with a smaller central value, 
corresponding to a |a downward fluctuation of the '98 measurement. For both hypothetical 
data sets, the global CL is greater by a factor two than the CL of the current data. 

To understand how correlates with the asymmetry measurements we also exhibit 
the corresponding fits in which either A''^^ or Alh are excluded. In the current data, the CL 
increases appreciably, to 0.51, if Ap^ is excluded but much less if Alr is excluded, refiecting 

^ The degrees of freedom decrease from 14 to 12 because we follow the recent practice of the EWWG||l| in 
consolidating the LEP II and FNAL measurements into a single mw measurement and the two t polarization 
measurements into a single quantity that we denote A^^t- The same set of measurements is tracked for all 
three years. 
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the larger pull of Ap^ in the SM fit. In the two hypothetical scenarios the CL increases 
comparably whether Alr or Ap^ is excluded. 

There are two conclusions from this exercise. First, the evolution of the mw measure- 
ment contributes as much to the marginalization of the global fit as does the evolution of 
the asymmetry measurements. Second, at its current value and precision mw tilts the SM 
fit toward Alr and xjy[74i:], while tagging Ap^ and x^yf^lij] as 'anomalous'. The reason for 
this "alliance" of mw and x^yfA^] will become clear in section 5, where we will see that mw 
and x^f^li] favor very hght values of the Higgs boson mass, far below the 114 GeV lower 
limit, while favors much heavier values, far above 114 GeV. 

Returning to table 2.2, we also see the effect on the global fit of the new result from 
NuTeV. In the '98 and '01 data sets, NuTeV had little effect on the global CL. In the 
current data set, because of its increased precision and central value, it causes the CL to 
decrease from an already marginal 0.10 to a poor 0.010. The low confidence level of the global 
SM fit is then due in roughly equal parts to 1) the discrepancy between the x\Y[Ai^]^mw 
aUiance versus ^'^[A/^], and 2) the NuTeV result. We will refer to and N] 

as "anomalous" simply as a shorthand indication of their deviation from the SM fit, with no 
judgement intended as to their bona fides. 
3. Interpreting the Discrepancies 

In this section we wish to set the context for the fits to follow by briefly discussing 
the three generic explanations of the discrepancies in the SM fit reviewed in the previous 
section. They are statistical fluctuation, new physics, and underestimated systematic error. 
Combinations of the three generic options are also possible. 

3.1 Statistical Fluctuations 

One or both anomahes could be the result of statistical fluctuations. However, if the data 
is to be consistent not only with the global fit but also with the lower limit on mH from the 
Higgs boson searches as discussed in sections 5 and 6, it is necessary that both anomalous 
and non-anomalous measurements have fluctuated. If only the anomalous measurements 
were to have fluctuated, the global flt would improve but the conflict with the lower limit 
on mn would be exacerbated. 

A high energy physics sage is reputed to have said, only partly in jest, that "The 
confidence level for 3a is fifty-fifty." The wisdom of the remark has its basis in two different 
phenomena. First, at a rate above chance expectation, many unusual results are ultimately 
understood to result from systematic error — this possibility is discussed below and its 
implications are explored in sections 5 and 6. Second, estimates of statistical significance are 
sometimes not appropriately defined. For instance, when a 3a "glueball" signal is discovered 
over an appreciable background in a mass histogram with 100 bins, the chance likelihood is 
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not the nominal 0.0027 associated with a 3a fluctuation but rather the complement of the 
probability that none of the 100 bins contain such a signal, which is 1 - 0.9973^°° = 0.24. 
The smaller likelihood is relevant only if we have an a priori reason to expect that the signal 
would appear in the very bin in which it was discovered. 

In assessing the possibihty of statistical fluctuations as the explanation of the poor SM 
flt, it should be clear that the global flt CL's are appropriately deflned, reflecting statistical 
ensembles that correspond to replaying the previous decade of experiments many times over. 
In particular, the CL's of the global flts are hke the glueball example with the signiflcance 
normalized to the probability that the signal might emerge in any of the 100 bins, as shown 
explicitly below. 

Table 3.1 summarizes fi^s of four different data sets, in which none, one, or both 
sets of anomalous measurements are excluded. Consider for instance flt B in which only 
x^[ v N\ is excluded, with CL = 0.10. In that flt, consisting of 16 measurements, the only 
signiflcantly anomalous measurement is A^fbi with a pull of 2.77, for which the nominal CL 
is 0.0056. The likelihood that at least one of 16 measurements will differ from the flt by 
> 2.77(7 is then 1 -0.9944^^ = 0.09, which matches nicely with the CL of 0.10. Similarly, 
in flt C which retains NuTeV while excluding the hadronic asymmetry measurements, the 
outstanding anomaly is v 'N\ with a pull of 3.0, and the probability for at least one 
such deviant is 1 - 0.9973^"^ = 0.04, compared to the ^ CL = 0.05. Finally, for the full data 
set, flt A, the outstanding anomalies are x^f^^i^^ N\ with a pull of 3.0 and with a pull of 
2.55. In that case we ask for the probability of at least one measurement diverging by > 3. Oct 
and a second by > 2.55ct, which is given by 1 - 0.9973^^ - 17(1 - 0.9973) (0.9892)^6 = 0.006, 
compared to the x^ CL of 0.010. 

We see then that the x^ CL's appropriately reflect "the number of bins in the his- 
togram," and that the poor CL's of these SM flts are well accounted for by the appropri- 
ately deflned probabilities that the outlying anomalous measurements could have occurred 
by chance. The nominal confidence levels of the global flts are then reasonable estimates 
of the probability that statistical fluctuations can explain the anomalies, which we may 
characterize as unlikely but not impossible. Only fit D, with both anomalies removed, has 
a robust confidence level, CL = 0.65. We refer to fit D as the "Minimal Data Set." The 
results and pulls for this fit arc shown in table 3.2. 

This discussion does not refiect the fact that the anomalous measurements are all within 
the subset of measurements that dominate the determination of mu- In that smaller subset 
of measurements, the significance of the anomalies is not fully refiected by the global CL's. 
As concerns the reliability of the fits of mj^, there is a clear a priori reason to focus on the 
mj^-sensitive measurements. We therefore also consider fits in which the observables Oi in 
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equation (4.1) below are restricted to the measurements which dominate the determination 
of rriH- These are the six asymmetries, m\Y, ^z, Ri, and A^]. (The m/^-insensitive 

measurements which are omitted from these fits are ah, Rb, Rc, Ab, and Ac.) The resuhs 
of the corresponding fits, A' - D', are tabulated at the bottom of table 3.1. Except for 
the minimal data sets, D and D', in every other case the fit restricted to m^f-sensitive 
measurements has an appreciably smaller confidence level than the corresponding global 
fit. In addition to the problems of the global fits, the poorer consistency of this sector 
of measurements provides another cause for concern in assessing the reliability of the SM 
prediction of mj^. 

3.2 New Physics 

Each anomaly could certainly be the result of new physics. The NuTeV experiment 
opens a very different window on new physics than the study of on-shell Z boson decays at 
LEP I and SLC.[|ri[| For example, a Z' boson mixed very little or not at all with the Z boson, 
could have little effect on on-shell Z decay but a big effect on the NuTeV measurement, which 
probes a space-like region of four- momenta centered around ~ —20 GeV^. The strongest 
bounds on this possibility would come from other off-shell probes, such as atomic parity 
violation, e+e~ annihilation above the Z pole, and high energy pp collisions. 

New physics could also affect the hadronic asymmetry measurements. Here we can imag- 
ine two scenarios, depending on how seriously we take the clustering of the three hadronic 
asymmetry measurements. Taking it seriously, we would be led to consider leptophobic Z' 
models [|I^], as were invoked to explain the Ri, anomaly, which was subsequently found to 
have a systematic, experimental explanation. 

Or we might regard the clustering of the three hadronic measurements as accidental. 
Then we would be led to focus on A^p^, by far the most precise of the three hadronic 
asymmetry measurements, and we could arrive at acceptable global fits by assuming new 
physics coupled predominantly to the third generation quarks. New physics would then 
account for the anomaly, with an additional effect on the less precisely measured jet 
charge asymmetry, Qfb. The third generation is a plausible venue for new physics, since the 
large top quark mass suggests a special connection of the third generation to new physics 
associated with the symmetry breaking sector. 

Since Rb oc g^L + qIr while A^p^ oc glj^ — glj^, and because Rb, which is more precisely 
measured than A^^, is only ~ la from its SM value, some tuning of the shifts 5gbL and 5gbR 
is required to fit both measurements. The right-handed coupling must shift by a very large 
amount, with 5gbR ^ 5gbL and 6gbR^0.1gbR. An effect of this size suggests new tree-level 
physics or radiative corrections involving a strong interaction. 

Examples of tree-level physics are Z — Z' mixing or b — Q mixing. A recent proposal 
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to explain the anomaly embeds a Z' boson in a right-handed SU{2) extension of the 
SM gauge group in which the third generation fermions carry different SU{2)ji charges than 
the first and second generations. Z' bosons coupled preferentially to the third generation 
are generic in the context of topcolor models. [|1^] An explanation hj b — Q mixing, requires 
Q to be a charge —1/3 quark with non-SM weak quantum numbers; this possibility has been 
explored in the context of the latest data in reference and previously in 

If new physics explains the A^p^ anomaly, it must also affect Af,. If we use eq. (2.3) 
with the factor Af, = 0.1501(17) taken from the three leptonic asymmetry measurements 
(assuming lepton universality), we find that the experimental value A'^p^ = 0.0994(17) implies 
AblA^ps] = 0.883(18), which is 2.89a from ^[SM] = 0.935, CL = 0.004. However Ab is 
measured more directly at SLC by means of A'^^^pj^, the front-back left-right asymmetry. 
In the Summer of 1998 that measurement yielded Ab[Apppp] = 0.867(35), lower by 1.9cr 
than Afe[SM], lending support to the new physics hypothesis. But the current measurement, 
M^fblr] = 0.922(20), is only 0.6a below ^^[SM]. It no longer bolsters the new physics 
hypothesis but it is also not grossly inconsistent with ^^[A^^], from which it differs by lAAa, 
CL = 0.15. Combining Ab[A^pp] and ^[A^BLij] we find AbiA^pppj^ © ^fbla] = 0.900(13), 
which differs from the SM by 2.69cr, CL = 0.007. Thus while the A^pppj^ measurement no 
longer supports the new physics hypothesis, it also does not definitively exclude it. 

If either anomaly is the result of new physics, the SM fails and we cannot predict the 
Higgs boson mass until the nature of the new physics is understood. New physics affecting the 
NuTeV measurement and/or the hadronic asymmetry measurements will certainly change 
the relationship between those observables and the value of rriH, and could affect other 
observables in ways that change their relationships with mn- 
3.3 Systematic Uncertainties 

The two 3(7 anomalies each involve subtle systematic issues, having to do both with 
performance of the measurements and the interpretation of the results. With respect to 
interpretation, both use nonperturbative QCD models with uncertainties that are difficult 
to quantify. In both cases the experimental groups have put great effort into understanding 
and estimating the systematic uncertainties. Here we only summarize the main points. 

The central value for the NuTeV SM result is x^i^'u^ N) = 0.2277±0.0013(statistical) ± 
0.0009(systematic). The estimated systematic error consists in equal parts of an experimen- 
tal component, ±0.00063, and a modelling component, ±0.00064. Uncertainty from the z/g 
and Vf. fluxes makes the largest single contribution to the experimental component, ±0.00039, 
with the remainder comprised of various detector-related uncertainties. The modelling un- 
certainty is dominated by nonperturbative nucleon structure, with the biggest component, 
±0.00047, due to the charm production cross section. 
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Two possible nonperturbative effects have been considered (see ^7\, 0, |T8[ and references 
therein). One is an asymmetry in the nucleon strange quark sea, / x{s{x) — s{x)) dx 0. 
Using dimuons from the separate u and 17 beams, the NuTeV collaboration finds evidence 
for a —10% asymmetry within the NuTeV cross section model. If truly present, it would 
increase the discrepancy from 3.0 to 3.7cr.[|l^ For consistency with the SM, an asymmetry 
of ^ +30% would be needed. 

A second possible nonperturbative effect is isospin symmetry breaking in the nucleon 
wave function, dP{x) ^ Studies are needed to determine if structure functions can 

be constructed that explain the NuTeV anomaly in this way while maintaining consistency 
with all other constraints. A negative result could rule out this explanation, while a positive 
result would admit it as a possibility. Confirmation would then require additional evidence. 

The 3cr discrepancy between x^yfA^] and also raises the possibility of subtle 

systematic uncertainties. The determination of xj^fAi] from the three leptonic measure- 
ments, Ai^ij, v4^^, and Ag.r, involve three quite different techniques so that large, common 
systematic errors are very unlikely. The focus instead is on the hadronic measurements, 
J^p^, Ap^, and Qfb- In these measurements, b and c quarks are mutual backgrounds for 
one another. The signs of both the A''p^ and Ap^ anomalies are consistent with misidenti- 
fying b ^ c, although the estimated magnitude] 19] of the effect is far smaller than what is 
needed. QCD models of charge flow and gluon radiation are a potential source of common 
systematic uncertainty for all three measurements. The two heavy flavor asymmetries, A^p^ 
and Ap^, have the largest error correlation of the heavy flavor Z-pole measurements, quoted 
as 16% in the most recent analysis. |l|] 

Since x'^^fAj:/] is dominated by the greatest concern is the systematic un- 

certainties of A^pQ. The combined result of the four LEP experiments is A^p^ = 0.0994 ib 
0.00157(statistics) ± 0.00071 (systematic). The systematic component arises from an "inter- 
nal" (experimental) component of ±0.00060 and a "common" component of ±0.00039, where 



the latter is dominated by the ±0.00030 uncertainty ascribed to QCD corrections. [[19 
3.4 Summary 

It is not now possible to choose among the three generic explanations of the anomalies, 
except to say that statistical fluctuations are unlikely per the nominal CL of the global fit. 
Bigger systematic errors could rescue the SM fit but would have to be much bigger than 
current estimates. Rather than further refinement of existing error budgets, this probably 
means discovering new, previously unconsidered sources of error. In this paper we focus on 
the systematic error hypothesis, not because we think it the most likely explanation — we 
do not — but, assuming it to be true, to see if it can really reconcile the SM with the data. 
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4. Methods 

In this section we describe the methodology of the SM fits. We also discuss the choice 
of observables, which differs slightly from the choice made in . 

We use ZFITTER v6.30p0[ to compute the SM electroweak radiative corrections, with 
results that agree precisely (to 2 parts in 10^ or better) with those obtained in The input 
parameters are m^, rrit, the hadronic contribution to the renormalization of a at the Z-pole, 
Aa^^mz), the strong coupling constant at the Z-pole, as{mz), and the Higgs boson mass, 
tuh- For any point in this five dimensional space ZFITTER provides the corresponding SM 
values of the other observables, Oi, listed in table 2.1. 

To generate the distributions we scan over the four parameters, rrit, Aa5{mz), 
asimz), and mn- For a specified collection of observables Oi, we then have 

The experimental values are given in table 2.1. 

We do not scan over mz but simply fix it at its central experimental value. Because 
mz is so much more precise than the other observables, it would contribute negligibly to 
if we did scan on it. We have verified this directly by performing fits in which it was varied, 
with only negligible differences from the fits in which it is fixed at the central value. This 
can also be seen in the global fits reported by the EWWG, in which the pull from mz is 
invariably much less than 1. In this case, inclusion of mz has no effect on the CL of the 
fit, because (1) the contribution to is negligible and (2) the scan on mz has no effect on 
the number of degrees of freedom since it is both varied and constrained. Since it has little 
effect, we choose not to scan on mz in order to facilitate the numerical calculations. 

For Aa^{mz) we use the determination of which incorporates the most recent e+e^ 
annihilation data and is also the default choice of |]l[. In |[T^ we also presented results for 
four other determinations of Aa^{mz)- 

For the global fits asijnz) is left unconstrained, as is also done in 0, because the Z-pole 
SM fit is itself the most precise determination of ^^(m^). For the fits which consider more 
limited sets of observables, we use the following rule: if at least two of the three observables 
which dominate the the determination of ^^(m^) (these are T z-, Ri, and ah) are included in 
the fit, is unconstrained as in the global fits. Otherwise we constrain it to 0.118(3). In 
any case, because the m//-sensitive observables are predominantly ag-insensitive, the results 
we obtain for mn depend very little on the details of how as is specified. 

The fits also include the important correlations from the error matrices presented in 
tables 2.3 and 5.3 of [^. We retain the correlations that are > 0.05 in the 6x6 heavy fiavor 
error matrix, for Ap^, Ap^, Aj,, A^ Rb, and Re- Similarly we retain correlations > 0.05 in 
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the 4x4 correlation matrix for Tz, (Jhi Ri, ci-nd A^p^. These correlations shift the value of 
by as much as 2 units. 

Our global fits differ slightly from the all-data fit of fl]], principally because we use the 
set of measurements the Electroweak Working Group (EWWG) used through Spring 2001 
but not two measurements that were subsequently added by the EWWG. From Summer 2001 
the EWWG all-data fit included the Cesium atomic parity violation (APV) measurement, 
and in Spring 2002, the W boson width, Tw, was included. 

Reference uses a 2001 determination of the Cesium APV measurement that has 
recently been superseded by newer results from the same authors|2^. With the average 
value from the more recent study, Qw = —72.18(46) (with experimental and theoretical 
errors combined in quadrature), and the SM value from [|l|, the pull is 1.6, rather than 0.6 
as quoted in |1[|. The effect on the all-data global fit in is to change x^/^ = 28.8/15, 
CL = 0.017 to x^/N = 30.8/15, CL = 0.009, decreasing the confidence level by a factor 
two. Rather than use the updated value, we choose not to include the Cs APV measurement 
since the theoretical systematic uncertainties are still in flux. 

We choose not to include Tw because it has not yet attained a level of precision precision 
approaching that of the other measurements in the fit. For instance, Tz is 30 times more 
precise than Tw, so that Tw has 1/900'th the weight of Tz in the global fits. At the current 
precision it has no sensitivity to new physics signals of the order of magnitude probed by 
the other observables in the fit. Its effect on the fits of mn, which are the principal concern 
of this work, is completely negligible. 

In any case, the decision not to include the and APV measurements does not have 
a major effect on our results. In particular, the effect on the Higgs boson mass predictions 
is negligible. Furthermore, our all-data fit, with CL = 0.010, has a very similar confidence 
level to that of the all-data fit of ||l|, CL = 0.009, if the APV determination is updated to 
reflect 

We closely reproduce the results of when we use the same set of observables. For 
instance, adding F^y and Qw{Cs) as specified in |I| to the observables in our global fit, table 
2.1, we obtain x^/N = 28.7/15, CL = 0.018, compared to x^/N = 28.8/15, CL = 0.017 
from the corresponding (all-data) fit of [0]. 

5. The x^ Fits 

In this section we present several SM fits of the data, using to estimate global CL's 
and to obtain the constraints on the Higgs boson mass. The Ax^ method, used also 
in |l|, is defined as follows. Let ttt-min be the value of mn at the minimum, and let niL 
be an arbitrary mass such that rriL > mMiN- Then the confidence level CL{mH > rni) is 
one half of the confidence level corresponding to a distribution for one degree of freedom. 
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with the value of the ^ distribution given by Ax^ = X^("^l) ~ X^("^min)- We consider both 
global fits and fits restricted to the mi^-sensitive observables. 

In addition to mt and Acts, which are input parameters to the ZFITTER calculations, 
the observables with the greatest sensitivity to mu are from the six asymmetry measure- 
ments, mwi ^z, Ri-i and x^{^i^ N]. It is useful to consider the domains in mn favored by 
these observables, in order to understand the "alliances" (see section 2) that shape the global 
fit, and to understand the consistency of the fits with the LEP II lower limit on mu- In 
fits restricted to mij-sensitive obervables, is given by eq. 4.1, where the Oi are restricted 
to the m//-sensitive observables under consideration. In addition, for fits containing fewer 
than two of the three as-sensitive observables — F^, i?;, and Uh — we constrain as{mz) by 
including it with the Oi in equation (4.1) as discussed in section 4. 

The experimental quantity that currently has the greatest sensitivity to ruH is x\^, 
determined from the six asymmetry measurements. Figure 1 displays the distributions of 
the three individual leptonic measurements and the combined result from all three, x^;^^[74i]; 
it shows that the upper limit is dominated by Alr. The central value, ruH — 55 GeV, 
symmetric 90% confidence interval, and hkelihood CL{mH > 114) are given in table 5.1. 
Note that the 95% upper hmit is ttih < 143 GeV. 

Figure 2 shows that the x^ distribution from x^;^^[74jj] is completely dominated by the b 
quark asymmetry, Ap^. The central value is mn — 410 GeV and the 95% lower hmit is 145 
GeV, just above the 95% upper hmit from The 95% upper limit extends above 1 

TeV. Figure 3 shows the x^ distributions of both x'^yfA^] and a;(y[/lj:/], with the respective 
symmetric 90% confidence intervals indicated by the dot-dashed horizontal lines. 

It is also interesting to isolate the effect of the W boson mass measurement, because 
it is the second most important quantity for fixing niH, and because it is a dramatically 
different measurement with a completely different set of systematic uncertainties. Figure 4 
shows the distribution from myy. The central value is run = 23 GeV and the 95% upper 
limit is 121 GeV. Both mw and xjyf/l^,] then favor very light values of mn- This is the 
basis of the "alliance" between ^'^[Al] and mw, discussed in section 2, that pushes Ap^ to 
outlyer status and contribTites to the marginal confidence level of the SM fit. 

The two other Z-pole observables with sensitivity to mn are Fz and Ri, which also 
involve different systematic uncertainties than the asymmetry measurements. The corre- 
sponding distributions are also plotted in figure 4, together with the combined distribu- 
tion for mw, and Ri. Ri also favors small mn, with its minimum off the chart below 
10 GeV. Tz is often represented as favoring mn — O(IOO) GeV, but we see in figure 4 and 
with the expanded scale in figure 5, that it actually has two nearly degenerate minima, at 



13 



about 16 and 130 GeV.j] In table 5.1 we see that the combined distribution of the three 
non-asymmetry measurements has a central value at 13 GeV and a 95% upper limit at 73 
GeV with CLimu) > 114) = 0.021 

Figure 6 shows the distribution from the NuTeV measurement, a;|^[^z/'' A^]. The 
minimum lies above 3 TeV and the 95% lower limit is at ~ 660 GeV. The SM fits the 
NuTeV anomaly by driving itlh to very large values, but the new physics that actually 
explains the effect, if it is genuine, would not be so simply tied to the symmetry breaking 
sector but might for instance reflect an extension of the gauge sector, with implications for 
the symmetry breaking sector that cannot be foretold. In any case, as discussed in section 
7, values of mn above ~ 1 TeV cannot be interpreted literally. 

It is striking that the measurements favoring mn in the region allowed by the direct 
searches are precisely the ones responsible for the large of the global fit. They favor values 
far above 100 - 200 GeV while the measurements consistent with the fit favor values far below, 
as shown in the bottom two lines of table 5.1. The fit based on x(y [A^] © mw ®Tz Q) Ri (fit 
D' in table 3.1) has < 106 GeV at 95% CL, while the fit based on x\^[Ah] ® x^fu^ N] 
has rriH > 220 GeV at 95% CL. The corresponding distributions are shown in figure 7. 

Next we consider the global fits that were discussed in section 3.1. The principal results 
are summarized in table 5.2. The "all-data" fit, fit A (shown in detail in table 2.1), closely 
resembles the all-data fit of |jl|, up to small differences arising from the slightly different 
choice of observables discussed in section 4. As summarized in table 5.2, the minimum 
is at mn = 94 GeV, with CL{x^) = 0.010. The 95% upper limit is mn < 193 GeV and 
the consistency with the search limit is CL{mH > 114 GeV) = 0.30. The distribution 
is shown in figure 8, where the vertical dashed line denotes the direct search limit and 
the horizontal dot-dashed line indicates the symmetric 90% CL interval. The combined 
likelihood for internal consistency of the fit and consistency between fit and search limit, 
defined in equation (1.1), is Pc = 0.0030. 

The "Minimal Data Set," fit D, with xjy [^z^'' A^] and omitted, is shown in detail 

in table 3.2, and the distribution is shown in figure 8. The minimum is at rriH = 43 GeV, 
with a robust confidence level CL(x^) = 0.65. But the 95% upper limit is ttih < 105 GeV 
and the consistency with the search limit is a poor CL{mH > 114 GeV) = 0.035. (The latter 
is nearly identical to the value 0.038 shown in table 5.1 for fit D', which is the corresponding 
fit restricted to m/^-sensitive observables, defined in section 3.1.) The combined likelihood 
for fit D is Pc = 0.023. 

In fit C with ^(^/[Aji/] omitted the minimum is at rriH = 45 GeV with CL{x^) = 0.066. 

wish to thank D. Bardin and G. Passarino, who kmdly verified this surprising feature, using, respec- 
tively, recent versions of ZFITTER[|o| and TOPAZO-H 
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The confidence level for consistency with the search limit is CL{mH > 114 GeV) = 0.047. 
The combined likelihood is Pc = 0.0031. 

Finally we consider fit B, with retained and x^fu^ N] omitted. Now the rriH 

prediction is raised appreciably, with the minimum at ttlh = 81 GeV, but the quality 
of the fit is marginal, with CL{)(^) = 0.10. The confidence level for consistency with the 
search limit is more robust, CL{mH > 114 GeV) = 0.26, and the combined likelihood is 
Pc = 0.026. 

The effect of the hadronic asymmetries and the NuTeV measurement is apparent from 
table 5.2. The NuTeV measurement diminishes CL(x^) by a factor 10, seen by comparing fit 
A with fit B and fit C with fit D, while its effect on CL{mH > 114) is modest. Consequently 
the NuTeV measurement also diminishes the combined likelihood Pc by an order of magni- 
tude. Comparing fit A with fit C or fit B with fit D, we see that the hadronic asymmetries 
also diminish CL(x^), by a factor ~ 7, but that they increase CL{mH > 114) by a nearly 
identical factor, so that they have little effect on Pc- 
6. "Bayesian" Maximum Likelihood Fits 

The Ax^ method for obtaining the confidence levels for different regions of rriH is poweful 
and convenient but, at least to this author, not completely transparent. We have therefore 
also approached these estimates by constructing likehhood distributions as a function of 
rriH, varying the parameters to find the point of maximum likelihood for each value of m^^. 
Assuming Gaussian statistics, the log likelihood is 

-logmmH))=Cx\mH), (6.1) 

so that the maximization of the likelihood is equivalent to the minimization of x^- The 
proportionality constant C is determined by the normalization condition for C 

The method is "Bayesian" in the sense that the domain of normalization and the measure 
are specified by a priori choices that are guided by the physics. The likelihood distribution 
is normalized in the interval 10 GeV < niu < 3TeV, 

f.mi/=3TeV 
/mH=10GeV 

The choice of measure is motivated by the fact that log(mi^) is approximately linearly 
proportional to the experimental parameters, which are assumed to be Gaussian distributed. 



/ d\og{mH)C{mH) = 1. (6.2) 

Jmw=10GeV 



such as nit, Aa^ and the various Oi — see for instance the interpolating formulae in 
The choice of interval is conservative in the sense that enlarging the domain above and below 
causes CL{mH > 114GeV) to be even smaller than the results given below. 

The normalized likelihood distributions for fits A and D, the all-data and Minimal Data 
Sets, are shown in figure 9, where we display both the differential and integrated distributions. 
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The confidence level CL{mii > 114) is the area under the differential distribution above 
114 GeV. For the Minimal Data Set the resuh is CL{mH > 114 GeV) = 0.030, in good 
agreement with the result 0.035 obtained from Ax^ in section 5. For the all-data set it is 
CL{mH > 114 GeV) = 0.25, compared with 0.30 from Ax^. 

It is clear from figure 9 that the likelihood distribution from the Minimal Data Set 
is vanishingly small at 3 TeV but has some support at 10 GeV. If we were to enlarge the 
domain in tuh both above and below, the effect would be to further decrease the likelihood 
for > 114 GeV. 

7. New Physics in the Oblique Approximation 

If we assume the Minimal Data Set, the contradiction with the LEP II lower limit on 
rriH is either a statistical fluctuation or a signal of new physics. Two recent papers provide 
examples of new physics that could do the job. Work by Altarelli et al.^E\ in the framework 



of the MSSM finds that the prediction for mn can be raised into the region allowed by the 
Higgs boson searches if there are light sneutrinos, ~ 55 - 80 GeV, light sleptons, ~ 100 
GeV, and moderately large tan P ~ 10. This places the sleptons just beyond the present 
experimental lower limit, where they could be discovered in Run II at the TeVatron. A 
second proposal, by Novikov et a/.,[0 finds that a fourth generation of quarks and leptons 



might also do the job, provided the neutrino has a mass just above mz/2. An illustrative 
set of parameters is niN — 50 GeV, uie — 100 GeV, mu + 171^ — 500 GeV, \mu — ttidI — 
GeV, and ^ 300 GeV. 

In this section we do not focus on any specific model of new physics but consider the 
class of models that can be represented in the oblique approximation [^, parameterized 



by the quantities S, T, U . |]^ The essential assumption is that the dominant effect of the new 
physics on the electroweak observables can be parameterized as effective contributions to the 
W and Z boson self energies. These contributions are not limited to loop corrections, since 
the oblique parameters can also represent tree level phenomena such a.s Z — Z' mixing. 
We will restrict ourselves to the S and T parameters, since they suffice to make the point 
that oblique new physics can remove the contradiction between the Minimal Data Set and 
the search limit, leaving mu as an essentially free parameter. We also show that S,T 
corrections do not improve the confidence levels of the global fits that include the anomalous 
measurements. 

For the observables Oj the oblique corrections are given by 

dQ, = Y.{A,S + B,T) (7.1) 



with Qi defined by Qi = Oi or Qi = In(Oi) as indicated in table 7.1 where Ai and Bi 
are tabulated. Since these are not SM fits, instead of N] the NuTeV experiment 
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is represented by the model independent fit to the effective couphngs g1 = g"^]^ + g^i and 
9%, — 9uR + 9dR^ which the experimental values from |^ are given in table 2.1. 

Figure 10 shows the S, T fit to the Minimal Data Set along with the SM fit with S = 
T = 0. The striking feature of the S,T fit is that is nearly fiat as a function of mn- 
There is therefore no problem reconciling the fit with the lower limit on m^^, and there is 
also no preference for any value of rriH- The fits are acceptable all the way to m^^ = 3 TeV, 
and the variation across the entire region is bounded by Ax^ < 1.2. Because the minimum 
is so shallow, it is not significant that it occurs at rriH = 17 GeV. The confidence level at 
the minimum is 0.51, which is comparable to the confidence level, CL = 0.65, of fit D, the 
corresponding SM fit. 

It is well known that arbitrarily large values of mn can be accomodated in S, T fits of 
the electroweak data. This can be understood as a consequence of the fact the SM fit of 
rriH is dominated by two observables, x\y and mw- Let m^fMIN], a;^y[MIN] and mvi/[MIN] 
be the values of itlh, x[y and mw at the minimum of the SM fit. The shifts and 
6mw induced in the SM fit by choosing a different value of mn m//[MIN], can then be 
compensated by choosing 5* and T to provide equal and opposite shifts, —6x\y and —6m\Y. 
Inverting the expressions from equation (3.13) of we have explicitly, in the approximation 
that we consider only and m\Y, 

S = --( Sx^w + 2x^J-^] (7.2) 

and 

r = --^ffeiv + ^). (T.3) 

a(l - Xw) \ mw ) 

For instance, for mn = 1 TeV the S, T fit to the Minimal Data Set (set D) shown in figure 
10 yields S,T = -0.22, +0.34 compared with S,T = -0.15, +0.22 from equations (7.2) and 
(7.3). The approximation correctly captures the trend though it differs by 30% from the 
results of the complete fit. The discrepancy refiects the importance of variations among 
parameters other than Xw and mw that are neglected in deriving (7.2) and (7.3). 

Values of mn above 1 TeV cannot be interpreted literally as applying to a simple 
Higgs scalar. For mn > 1 TeV symmetry breaking is dynamical, occurring by new strong 
interactions that cannot be analyzed perturbatively. [|] If the Higgs mechanism is correct, 
there are new quanta that form symmetry breaking vacuum condensates. Values of m^ 
above 1 TeV should be regarded only as a rough guide to the order of magnitude of the 
masses of the condensate-forming quanta. 

It is sometimes said that an SM Higgs scalar above ~ 600 GeV is excluded by the 



triviality bound, which is of order 1 TeV in leading, one loop order |PD[, refined to ~ 600 
GeV in lattice simulations. [0] The bound is based on requiring that the Landau singularity 
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in the Higgs boson self-coupling, A, occur at a scale ALandau that is at least twice the Higgs 
boson mass, ALandau '^Sm//, in order for the SM to have some minimal "head room" as an 
effective low energy theory. However, the conventional analysis does not include the effect on 
the running of A from the new physics which must exist at the Landau singularity. Although 
power suppressed, in the strong coupling regime, which is in fact the relevant one for the 
upper limit on uih, the power suppressed corrections can change the predicted upper limit 
appreciably, possibly by factors of order one.||32|| To take literally the 600 GeV upper bound 
from lattice simulations we in effect assume that the new physics is a space-time lattice. 
The bound cannot be known precisely without knowing something about the actual physics 



that replaces the singularity. The analysis in |32] is performed in the symmetric vacuum 
and should be reconsidered for the spontaneously broken case, but the conclusion is likely 
to be unchanged since it follows chiefly from the ultraviolet behavior of the effective theory 
which is insensitive to the phase of the vacuum. An SM scalar between 600 GeV and 1 TeV 
therefore remains a possibility. 

Figure 10 also displays the values of S and T corresponding to the minimum at each 
value of rriH- For itlh > 114 GeV the minima fall at moderately positive T and moderately 
negative S. Positive T occurs naturally in models that break custodial SU{2), for instance 
with nondegenerate quark or lepton isospin doublets. Negative 5* is less readily obtained but 
there is not a no-go theorem, and models of new physics with 5* < have been exhibited. |p3| 

We also consider a fit to the Minimal Data Set in which only T is varied with S held at 
S = 0. The result is shown in figure 11. The minimum falls at rriH = 55 GeV with CL = 0.56. 
The distribution at larger rriH is fiat though not as flat as the S,T fit. Moderately large, 
postive T is again preferred. From Ax^ we find that the confidence level for mn above the 
LEP II lower limit is sizeable, CL{mH > 114 GeV) = 0.21, and that the 95% upper limit is 
tuh < 460 GeV. 

Next we consider the S, T fit with the hadronic asymmetry measurements included, 
corresponding to SM fit B above. Shown in figure 12, it is not improved relative to the SM 
fit. The minimum is at mn = 15 GeV with x^/^ = 15/10 implying CL = 0.12. From 
Ax^ the probability for mn in the allowed region is a marginal CL{mH > 114 GeV) = 0.08, 
and the combined probability from equation (1.1) is Pc = 0.01. 

The all-data S, T fit, including both the hadronic asymmetry and NuTeV measurements, 
is shown in figure 13. In this case the S* = T = fit is not identical to SM fit A, since the 
NuTeV result is parameterized by gi and gn rather than Xyy[ v N\ as in the SM fit. The 
minimum of the S" = T = fit occurs at mu = 94 GeV, with x^ jN = 26/14 implying 
CL = 0.026. The S,T fit is actually of poorer quality: the shallow minimum is at tuh = 29 
GeV with x'^/N = 25.7/12 implying CL = 0.012. 
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8. Discussion 

Taken together the precision electroweak data and the direct searches for the Higgs 
boson create a complex puzzle with many possible outcomes. An overview is given in the 
"electroweak schematic diagram," figure 14. The diagram illustrates how various hypotheses 
about the two 3a anomalies lead to new physics or to the conventional SM fit. The principal 
conclusion of this paper is refiected in the fact that the only lines leading into the 'SM' 
box are labeled 'Statistical Fluctuation.' That is, systematic error alone cannot save the 
SM fit, since it implies the confiict with the search limit, indicated by the box labeled 
CL{mH > 114) = 0.035, which in turn either implies new physics or itself refiects statistical 
fiuctuation. This is a consequence of the fact that the combined probability Pc defined in 
equation (1.1) is poor whether the anomalous measurements are included in the fit or not, 
as summarized in table 5.2. 

The 'New Physics' box in figure 14 is reached if either 3a anomaly is genuine or, con- 
versely, if neither is genuine and the resulting 96.5% CL confiict with the search limit is 
genuine. It is also possible to invoke statistical fiuctuation as the exit line from any of the 
three central boxes. However we have argued that the global confidence levels indicated 
for fits A and B are fair refiections of the probability that those fits are fiuctuations from 
the Standard Model. As such they do not favor the SM while they also do not exclude it 
definitively: "It is a part of probability that many improbable things will happen." 



The smoothest path to the SM might be the one which traverses the central box, fit 
B, and then exits via 'Statistical Fluctuation' to the SM. In this scenario nucleon structure 
effects might explain the NuTeV anomaly and the 10% confidence level of fit B could be 
a fiuctuation. This is a valid possibility, but two other problems indicated in the central 
box should also be considered in evaluating this scenario. First, the consistency of the mn- 
sensitive measurements is even more marginal, indicated by the 3.4% confidence level of 
fit B'. Second, the troubling 3a confiict {CL = 0.003) between the leptonic and hadronic 
asymmetry measurements is at the heart of the determination of mn- Thus even if we assume 
that the marginal CL of the global fit is due to statistical fiuctuation, the reliability of the 
prediction of tjih hangs on even less probable fiuctuations. As noted above, to be consistent 
with the search limit statistical fiuctuations must involve both the 'anomalous' hadronic 
asymmetry measurements and the measurements that conform to the SM fit, especially the 
leptonic asymmetry measurements and the W boson mass measurement. The confiict with 
the search limit would be greatly exacerbated if the true value of were equal to the 

present value of 

Since there are still some ongoing analyses of the hadronic asymmetry data, there may 
yet be changes in the final results, but unless major new systematic effects are uncovered 
the changes are not likely to be large. More precise measurements might be made eventually 
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at a second generation Z factory, such as the proposed Giga-Z project. However, to fully 
exploit the potential of such a facility it will be necessary to improve the present precision 
of Aa^lrriz) by a factor of ~ 5 or better, requiring a dedicated program to improve our 



knowledge of a{e^e^ HADRONS) below ~ 5 GeV. I^S] The W boson and top quark mass 
measurements will be improved at Run II of the TeVatron, at the LHC, and eventually at a 
linear e~^e~ collider. For instance, an upward shift of the top quark mass|^B[ or a downward 
shift of the W boson mass could diminish the inconsistency between the Minimal Data Set 
and the search limit, while shifts in the opposite directions would increase the conflict .0 

The issues raised by the current data set heighten the excitement of this moment in high 
energy physics. The end of the decade of precision electroweak measurements leaves us with 
a great puzzle, that puts into question the mass scale at which the physics of electroweak 
symmetry breaking will be found. The solution of the puzzle could emerge in Run II at the 
TeVatron. If it is not found there it is very likely to emerge at the LHC, which at its design 
luminosity will be able to search for the new quanta of the symmetry breaking sector over 
the full range allowed by unitarity. 
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Table 2.1. SM All-data fit (fit A). Experimental values for the model-independent parameters 
gliJ^v N] and g^iv N] are given for completeness but are not used in the SM fit. 





Exp er imcnt 
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0.2277 (16) 


0.2227 


3.0 


rut 


174.3 (5.1) 


175.3 


-0.2 


Att5(m|) 


0.02761 (36) 


0.02768 


0.2 


as{mz) 




0.1186 




ruH 




94 




9lt^^ N] 
9lt^^ N] 


0.3005 (14) 






0.0310 (11) 







24 



Table 2.2. Evolution of the electroweak data. As noted in the text, the same data is tracked 
for the three data sets though, following |jl|, it is grouped into fewer degrees of freedom in 
the Spring '02 data set. 



Summer '98 Spring '01 Spring '02 



x\y [Ah] 
CL{Al © Ah) 

mw 

X^N (no x^^fu N] 
CLixVN] 



0.23128 (22) 
0.23222 (33) 

0.02 

0.25 
80.410 (90) 

13.8/14 
0.46 



0.23114 (20) 
0.23240 (29) 
0.0003 
0.02 
80.448 (34) 

24.6/14 
0.04 



0.23113 (21) 
0.23220 (29) 
0.003 
0.06 
80.451 (33) 

18.4/12 
0.10 



X 



Pu\\{x^[^u N]) 
CLlxVN] 



0.2254 (21) 

1.1 

15/15 
0.45 



0.2255 (21) 

1.2 
26/15 
0.04 



0.2277 (16) 

3.0 
27.7/13 
0.01 



Table 2.3. "What if?": role of mw in shaping the global fit. The first column reflects actual 
current data with [ A^] omitted. In the second and third columns mw is assigned 
hypothetical values as described in the text, while other measurements are held at their 
Spring '02 values. In each case the effect of omitting Ap^ or Alr is also shown. 





mw['02] 


mw['98] 


AmvK['02] 




80.451 (33) 


80.410 (90) 


80.370 (33) 


XV12, CL 


18.4, 0.10 


15.2, 0.23 


15.3, 0.23 










xVll, CL 


10.2, 0.51 


9.0, 0.62 


9.8, 0.55 


OR 








-Alr 








xVll, CL 


15.7, 0.15 


10.2, 0.51 


10.0, 0.53 
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Table 3.1. Results for global fits A - D and for the corresponding fits restricted to tjih- 
sensitive observables, A' - D'. 





All 


-x^tv N] 


All 


A 

XV = 27.7/13, CL = 0.010 


B 

18.4/12, 0.10 


-x\^[Ah\ 


C 

17.4/10, 0.066 


D 

6.8/9, 0.65 


mif-sensitive only: 
All 


A' 

24.3/8, 0.0020 


B' 

15.2/7, 0.034 




C 

13.8/5, 0.017 


D' 

3.45/4, 0.49 



Table 3.2. SM fit D, to Minimal Data Set, with Xy^[ v N\ and three hadronic asymmetry 
measurements excluded. 





Experiment 


SM Fit 
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Table 5.1. Predictions for uih from various restricted sets of mij-sensitive observables. The 
value of ruH at the minimum is shown along with the symmetric 90% confidence interval 
and the likehhood for m^^ > 114 GeV. Values indicated as 10— or 3000+ fall below or above 
the interval 10 < m^j < 3000 GeV within which the fits are performed. 





rriH (GeV) 90% CL CL{mH > 114) 


mw 

mw ®Tz®Ri 

x^i^^ N] 


55 16<mH< 143 0.10 
410 145 <mH < 1230 0.98 
23 10- <mH< 121 0.059 
13 10- < < 73 0.021 

3000+ 660 <mH < 3000+ 0.996 


x\v[AL]®mw®Tz®Ri 


37 11 < niH < 106 0.038 
600 220 <mH< 1690 0.995 



Table 5.2. Confidence levels and Higgs boson mass predictions for global fits A - D. Each 
entry shows the value of niH at the minimum, the symmetric 90% confidence interval, 
the confidence level, the confidence level for consistency with the search hmit, and the 
combined likelihood Pc from equation (1.1). 
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CL(mH > 114) = 0.30 
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CL(x^) = 0.65 




CL{mH > 114) = 0.047 


CL{mH > 114) = 0.035 




Pc = 0.0031 


Pc = 0.023 
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Table 7.1 Coefficients for the oblique corrections as defined in equation (7.1). 
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0.000659 


-0.000468 


Ac 


-0.0125 


0.00886 


At, 


-0.00229 


0.00163 


9l 


-0.00268 


0.00654 


9l 


0.000926 


-0.000198 



28 



Leptonic Asymmetries 



X 
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mH(GeV) 



Figure 1: distributions as a function of itlh from the leptonic asymmetry measurements. 
The dashed, dot-dashed, and dotted hues are obtained from Alji, A^p^, and A^^t- respectively. 
The solid line is the combined fit to the three asymmetries. 
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Hadronic Asymmetries 




mH(GeV) 

Figure 2: distributions from the hadronic asymmetry measurements. The dashed, dot- 
dashed, and dotted Unes are obtained from A^p^, Ap^, and Qfb respectively. The sohd hne 
is the combined fit to the three asymmetries. 
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Leptonic &c Hadronic 




10 20 50 100 200 500 1000 2000 



mH(GeV) 

Figure 3: distributions from the leptonic (left) and hadronic (right) asymmetry measure- 
ments. The dot-dashed lines indicate the respective symmetric 90% CL intervals. 
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Figure 4: distributions from non-asymmetry measurements. The dashed, dot-dashed, and 
dotted hnes are obtained from mw, ^z, and Ri respectively. The sohd hne is the combined 
fit to the three measurementss. The dot-dashed hnes mark the 95% CL upper hmits for the 
combined distribution and for the distribution obtained from mw alone. 
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Figure 5: distribution from V z with expanded scale. 
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NuTeV 

n I ' ' I ""I n 
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Figure 6: distribution from x^fu N]. 
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mTT — Sensitive Observables 
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m„(GeV) 

Figure 7: distributions from the m/^-sensitive observables. The distribution on the left is 
a fit to the D' data set, i.e., restricted to the m//-sensitive observables Alr, A^p^, Ae^r, ^w, 
T Zi and i?/. The distribution on the right is a fit to the remaining mj:/-sensitive observables: 
A^p^, Ap^, Qfb, and a;j^[*^t^'' N]. The dot-dashed lines indicate the respective symmetric 
90% CL intervals. 
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All-Data (A) &c Minimal Set (D) 




10 20 30 50 70 100 200 



mH(GeV) 

Figure 8: Ax^ for two SM global fits. The Minimal Data Set, fit D, is on the left and the 
all-data set, fit A, is on the right. The 90% symmetric confidence intervals are indicated by 
the horizontal dot-dashed hne. The vertical dashed line denotes the 95% CL lower hmit on 
niH from the direct searches. 
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All Data (A) & Minimal Set (D) 




mH(GeV) 

Figure 9: Differential and integrated likelihood distributions for the Minimal Data Set (set 
D, solid hues) and the all-data set (set A, dotted hues). The vertical dashed line denotes 
the 95% CL lower limit onmn from the direct searches, and the horizontal dot-dashed hne 
indicates the 5% hkelihood level. 
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Minimal Data Set 



C\2 




10 20 50 100 200 500 1000 2000 



mH(GeV) 

Figure 10: distributions (solid lines) for the SM and S,T fits to the Minimal Data Set 
(data set D). The corresponding values of S (dashed line) and T (dot-dashed line) are read 
to the scale on the right axis. 
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Figure 11: distributions (solid lines) for the SM and T-only fits to the Minimal Data Set 
(data set D). T (dot-dashed line) is read to the scale on the right axis. 
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Fit B Data Set 



C\2 
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mH(GeV) 

Figure 12: distributions (solid lines) for the SM and S,T fits to data set B, i.e., including 

the hadronic asymmetry measurements but not a;^[''t^^ N]. S (dashed line) and T (dot- 
dashed line) are read to the scale on the right axis. 
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All— data 



CV2 




10 20 50 100 200 500 1000 2000 



mH(GeV) 

Figure 13: distributions (solid lines) for S,T fits to the all-data set, data set A. The 
distribution for 5" = T = is not equivalent to the SM fit since it uses the model-independent 
NuTeV fit to gL,R as discussed in the text. S (dashed hne) and T (dot-dashed hne) are read 
to the scale on the right axis. 
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Figure 14: Electroweak schematic diagram. 
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